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Global MRI with Braginskii viscosity in a galactic profile 
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ABSTRACT 

We present a giobal-in-radius linear analysis of the axisynimetric magnetorotational 
instability (MRI) in a coUisional magnetized plasma with Braginskii viscosity. For a 
galactic angular velocity profile 57 we obtain analytic solutions for three magnetic field 
orientations: purely azimuthal, purely vertical and slightly pitched (almost azimuthal). 
In the first two cases the Braginskii viscosity damps otherwise neutrally stable modes, 
and reduces the growth rate of the MRI respectively. In the final case the Bragin- 
skii viscosity makes the MRI up to 2^/2 times faster than its inviscid counterpart, 
even for asymptotically small pitch angles. We investigate the transition between the 
Lorentz-force-dominated and the Braginskii viscosity-dominated regimes in terms of 
a parameter Q.v^/ B'^ where is the viscous coefficient and B the Alfven speed. In 
the limit where the parameter is small and large respectively we recover the inviscid 
MRI and the magnetoviscous instability (MVI). We obtain asymptotic expressions for 
the approach to these limits, and find the Braginskii viscosity can magnify the effects 
of azimuthal hoop tension (the growth rate becomes complex) by over an order of 
magnitude. We discuss the relevance of our results to the local approximation, galax- 
ies and other magnetized astrophysical plasmas. Our results should prove useful for 
benchmarking codes in global geometries. 

Key words: instabilities - accretion, accretion discs - Galaxy: disc - MHD - magnetic 
fields - plasmas. 



1 INTRODUCTION 

In the formation of compact objects (stars, planets and 
black holes) from accretion discs, turbulence driven by the 
MRI, and possibly the MVI, offers a promising m echanism 
for the necessary angular momentum transport (IVelikhov 



tor tne necessary angular momentum transp oi 
19591 : IChandrasekhaiiil960l : iBalbus fc Hawlevl ll 



1991: Balbui 



200J). It has also been suggested that the observed velocity 



fluctuations ~ 6kms~ in parts of the interstellar medium 
(ISM) with low star formation rates may, in part, arise 
from this process (jSellwood fc Balbusll 19991 : iTamburro et al.l 
l2009l ). The evidence for this comes primarily from numerical 
simulations and a wide range of studies agree that weak mag- 
netic fields and outwardly decreasi ng angular velocity pro- 
files are an unstable combination (|Balbus fc Hawlevlll998l : 
lBalbus|[2003l ) 

The most illuminating explanation for this comes from 
a shearing sheet analysis in which the mean azimuthal flow 
of the differentially rotating disc is locally approximated by 
a cons tant angular velocity rotation p l us a linear velocity 
shear (|Goldreich fc Lynden-BeTil 1 19651 : lUmurhan fc RegevI 



E-mail: msr35@math.ucla.edu 



|2004 ). In the simplest possible setup, incompressible, 
isothermal, dissipationless, axisymmetric linear perturba- 
tions to a magnetic field with a weak vertical component, i.e. 
parallel to the rotation axis of the disc, are unstable when 
the angular velocity decreases away from the disc's centre. 
Azimuthal velocity perturbations to fluid elements at differ- 
ent heights, tethered to each other by the magnetic field, 
increase (decrease) their angular momentum. This causes 
them to move to larger (smaller) radii as dictated by the 
gravitational field which sets the mean flow. In the frame 
rotating at constant angular velocity, this motion deforms 
the tethering magnetic fleld, provided it is not too strong, 
and this induces a prograde (retrograde) Lorentz force on 
the outer (inner) element thus destabilising the system as 
it moves to yet larger (smaller) radii. This mechanism is at 
the heart of the MRI. 

However, although this model and description captures 
much of the essential physics, to fully understand the MRI, 
or at least the framework within which the shearing sheet 
should exist, a more nuanced approach is needed. In part, 
this is because the shearing sheet is formulated in a Carte- 
sian coordinate system where curvature terms that arise nat- 
urally from the cylindrical geometry of the accretion disc 
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are neglected. Indeed, in the local approximation, the over- 
stabilising effects of hoop tension (a curvature effect) associ- 
ated with the azimuthal magnetic field, are totally ignored. 
Furthermore, the model predicts that the fastest growing, 
and therefore most physically relevant, linear MRI modes 
have a homogeneous radial structure on the scale of the 
s hearing shee t in w hich the local approximation is made 
( iGuan et al.l (|2009l ) has shown the MRI is well localized 
in the non-linear regime). This means that the global disc 
structure, including boundary conditions, not captured by 
the local approximation may have a significant effect on 
these large scale modes in a way that c annot be deter- 
mined locally. Other lim itations exist too (|Knoblochl 1 19921 : 
iRegev fc Umurhanllioog ). 

The extent to which these limitations matter 
should, and under a variety of assumptions have, 
been investigated by global analyses that take into 
account the f ull radial structure of the disc and 
its b o undaries (iDubruUe fc KnoblochI 19931: Currv et al.l 



1994 ICurrv fc Pudritj 1 19951 : lOgilvie fc Pringld 1 19961 : 
Ogilvid 1998h, a nd specifically for the galaxy by 
Kitchatinov fc Riidigcr ( 2003 '). The conclusions of these in- 



vestigations largely confirms the local picture of a large ra- 
dial scale instability driven by the differential rotation of the 
disc. This suggests that whilst a local MRI analysis is gener- 
ally correct, its regime of validity must be checked globally. 
It is the purpose of this work to do just that for the MRI 
operating in a coUisional, magnetized (the ion cyclotron fre- 
quency oji 3> ion-ion collision frequency vu) plasma (like the 
ISM). 

In such a plasma iBraginskiil (jigesT l has shown that to 
lowest order in uji/vu, the deviatoric stress tensor is diag- 
onal and anisotropic. This leads to different parallel and 
perpendicular viscosities or, more fundamentally, pressures 
(and thermal conductivities]]) with respect to the local direc- 
tion of the magnetic field. Of the important physical conse- 
quences of this, it will be the effect of the Braginskii viscosity 
in the presence of a galactic shear flow that will concern us 
here. 

This is not a new topic and in recent years the study 
of magnetized accreti on discs has attracted attention in 
both the coUision less (|Quataert et al ] l2002l : ISharma et al.l 



20031. l2006l. l2007l) aii d coUisional reg ii nes (iBalbud |2Q04|: 
Islam fc BalbusI l2005l : iFerrard l2007i : [ Devlen fc Pekunlfil 



2010t). and a well developed code to simulate them now exists 



( Parrish fc Stondl2007l : Istone et al.|[2"008l ). However, a num- 
ber of fundamental questions remain unanswered. Primarily, 
what is the non-linear fate of the MRI in a coUisional mag- 
netized plasma? Does it transport angular momentum and 
if so, is the transporting stress primarily viscous. Maxwell 
or Reynoldffl? On what scales do the most unstable modes 
emerge and how does this vary with the system parame- 
ters? In what regime will a local analysis become untenable 



and global effects (either radia l or vertical) become impor- 
tant jCammie fc Balbudliggj )? What effect does the pres- 
ence, or absence, of a net vertical field have given Cowl- 
ing's anti-dynamo theorem and the diss i pative properti es 
of the Braginsku viscosity (|Moffatt|[T97i : lLvutikovl 120071 ')? 
Could viscous heating from the Braginskii viscosity lead 
to sec ondary magn e tized or unmagnetized thermal ins tabil- 
ities jBalbud [2OO1I : lQuataer'3 1200^ : [Kuuz et al.|[2oTll )? Do 
channel solutions , or so mething approaching them, emerge 
(|Goodman fc Xiil Il994l ')? If they do, the associated field 
growth will generate pressure anisotropies that could feed 
new parasitic instabilities such as the mirror. What would 
their effects be at this stage and in the inevitable turbu- 
lence where the mirror and fi rehose instabilities will both 
arise (|Schekochihin et al.ll2005h ? 

Addressing these questions will require a two-pronged 
approach involving both numerical and analytic studies. It 
may transpire that much of the existing work on the unmag- 
netized and coUisionless MRI is directly applicable, but this 
needs determining. As such we conduct a global linear sta- 
bility analysis for three separate background magnetic field 
orientations: purely azimuthal, purely vertical and pitched 
(magnetic field lines follow helical paths on cylinders of con- 
stant radius). We embed these in a galactic rotation profile. 
In agreement with earlier, local studies, we find that when 
the field has both a vertical and an azimuthal component, 
a linear instability wit h a real part up to 2y/2 times faste r 
than the MRI emerges (|Balbudl200i llslam fc Baibu3l2005l ) . 
In contrast to local studies we also find that it has a trav- 
elling wave component, and its growth rate depends on the 
viscous coefficient. 

In the presence of a vertical field, we also recover the 
standard inviscid MRI and show that upon introducing the 
Braginskii viscosity, its growth rate is always reduced. A 
similar effect is found for a purely azimuthal field where we 
find that the Braginskii viscosity damps modes that are, 
inviscidly, neutrally stable. 

The layout of this paper is as follows. In Section [21 we 
introduce and perturb a series of equilibrium solutions to 
the Braginksii-MHD equations and this forms the basis of 
our global stability analysis. Relegating the manipulation of 
the perturbed equations to Appendix [X] we obtain a single 
ODE governing the perturbed modes. We proceed by solv- 
ing this for azimuthal, vertical and pitched field orientations 
in Sections|3]|4]and[5]respectively. For each case we contrast 
the behaviour with and without Braginskii viscosity. In Sec- 
tion [S] we discuss the physical mechanism of the instability, 
where it may occur astrophysically, and the relation of our 
results to the local approximation. Finally, we conclude in 
Section [7] with some thoughts on open questions relating to 
magnetized astrophysical plasmas. 



^ In the presence of temperature gradients, a nisotropic the rmal 
conduction can lead to the magnetothermal j[Balbuj[2OO0l) 



and 



heat flux buoyancy instabilities jQuataerdbOoir 

Simulations in the coUisionless regime bv lSharma et al.l ||2006| ') 
shows that there is angular momentum transport and the 
anisotropic pressure constitutes a significant portion of the to- 
tal stress (~ Maxwell and 3> Reynolds). 



2 GLOBAL STABILITY ANALYSIS 
2.1 Governing equations 

The simplest set of equations that capture the physics of 
the coUisional magnetized MRI are those th ose of isother- 
mal i deal MHD with the Braginskii viscosity (|Lifshitz et al.l 
1 1984 ). Explicitly these equations are the momentum equa- 
tion, including the Braginskii stress tensor; the induction 
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equation that describes the evolution of the magnetic field; 
the incompressibility condition (because the perturbations 
are linear, the Mach number can always be made small 
enough to ensure this); and the stress tensor itself: 



du ^ 
— + u-\/u 
at 



-vn + -B • VB - V$D - V • T, 



dB 
~dt 



u VB = B ■ Vn. 



V ■ u = , 



T = (I - 36b) bb : Vu, 



(1) 

(2) 
(3) 

(4) 

where the constant mass density has been scaled out of the 
problem, u is the velocity field, B is the mass-density scaled 
magnetic field, i.e. the Alfven velocity, b = B/B is the direc- 
tion of the magnetic field, <&d is the gravitational potential, 
If = p + /2 is the total (gas plus magnetic) pressure, ':' 
is the full inner product, I the identity, and T is the full 
Braginskii stress tensor whose form we now explain. 

In a magnetized plasma the total pressure tensor P = 
pi -I- T (whose divergence appears in the momentum equa- 
tion) is given by 

P=pl + ^(l-3b6)(px-P||), 

where px,P|| are the perpendicular and parallel scalar pres- 
sures with respect to the magnetic field direction. When the 
plasma is also coUisional, the pressure anisotropy px— py can 
be related to the rat e-of-strain by a Chapman-Enskog style 
pertu rbation theory (|Brae:inskiill 19651 : IChapman fc Cowlind 
11973). 

Microphysically, the anisotropy arises from the conser- 
vation of the magnetic moment (first adiabatic invariant) of 
a gyrating particle in a magnetic field. However, collisions 
break this conservation and relax the anisotropy by pitch- 
angle scattering particles in velocity space (this dissipative 
process will turn out to be important in isolated field con- 
figurations). The competition between these two processes 
is governed by 



Vii{p± -P\\), 



d , , d In B 

-(px-P||)^3p^ 

where we have used the BGK operator to approximate the 
full collision operator. 

In the presence of time-varying magnetic fields, the 
pressure anisotropy tends to a steady state that tracks the 
fields' rate of change. It follows 



-P|| = 3^-3 - 



d\nB 



dt 



Si^Bbb : Vu, 



where vb ~ p/vu is the coefficient of the Braginskii viscosity, 
and we have used equation ((2} in the final equality. Equation 
Q follows directly. 

2.2 Equilibrium solutions 

Working in cylindrical polar coordinates {r,(j},z), we intro- 
duce equilibrium solutions to equations ([l))-© that describe 
a differentially rotating global shear flow constrained by 
gravity and threaded by a magnetic field that lies on cylin- 
ders of constant radius. Our equilibrium solutions are uni- 
form in z, and the plasma motion is restricted at an inner 
boundary of finite radius ro. (The validity of these assump- 
tions with be discussed in Section [6l2]) 



We allow gravity V'I'd to dictate the rotation profile 
of the equilibrium flow Q, — Q,o{r /ro)^'' where Q,o is the 
rotation frequency at the inner boundary and g is a dimen- 
sionless measure of the shear. Of interest to us here is the 
case of q — 1 that is both analytically treatable and phys- 
ically corresponds to a galactic disc where, unlike the Ke- 
plerian case of q — 3/2, the gravitational p otential of the 
dark matter halo sets the rotation proflle (|Rubin fc FordI 
Il97d : ISofue fc Rubiiil bOOlh . In modelling this we set the 
dark matter mass distribution (whose sole pu r pose is to ul- 
timately set the rotation profile) to a iMestej (| 19631 ) profile 
so, via Poisson's equation, V$d oc 1/r. In this case, the flow 
u — Q{r)re^ = ^Iqb^ does not vary with radius. 

We decompose the magnetic field into vertical and az- 
imuthal components, (so as to construct a time indepen- 
dent equilibrium, we neglect radial mag netic fieldil!!). That 
is, B — BzBz + B^e.^ = _B(sin^e, -f cosOe^) where 9 = 
arctan{Bz/Bif,) is the pitch angle of the magnetic field. We 
do not specify 9 yet, however, for mathematical simplicity, 
we demand that both B^ and B^ are independent of radius 
so 6 remains constant. Whilst this implies a vertical current 
oc 1/r is associated with B^, a simple super-galactic mag- 
netic field can account for B^. In all, our equilibrium fields 
take the form 



u = Sloe,* 



B — BzSz + B^e^, 



(5) 



and are constant in space and time. As is physically relevant, 
we restrict B /Q. < 1. 

It is a mathematically convenient feature of our equi- 
librium solutions that there is no evolution of the magnetic 
field strength. From equations @ and (O, T is absent from 
the unperturbed state and the system will be stable to pres- 
sure anisot ropy driven microscale in stabilities, e.g. firehose 
and mirror (|Schekochihin et al.l2005h . In contrast to the case 
where Laplacian viscosity (or indeed resistivity) is present, 
we can const ruct an ideal MHD solution independent of any 
radial flows (|Kersale et al.ll2004 ). 



2.3 Perturbed equations 

To determine the stability of this system we linearise equa- 
tions ([l|-(|4]) about ([5]) with axisymmetric velocity pertur- 
bations Su = 5u{r) exp[ifc2 -I- ^t] and similarly for the mag- 
netic and pressure fields. Here k is the wavenumber in the 
z direction and 7 the growth rate. We retain curvilinear 
terms from the cylindrical geometry but neglect self-gravity. 
We non-dimensionalise with respect to time-scales Slo and 
length-scales ro. 

As detailed in Appendix [X] the linearised equations 
combine into a single complex second order ordinary dif- 
ferential equation for 5ur, the Modified Bessel equation 



d 5ur 
dr-^ 



1 diSur 
r dr 



Sur = 0, 



(6) 



Although magnetic field configurations vary from galaxy to 
galaxy, they are commonly found tracing the spiral arms and 
therefore, in the plane of the disc, predominantly azimuthal. 
I Beck et al. I {mi) gives values for the mean in-plane field 
Br/B^ ~ 0.25 thereby justifying, to some degree, our neglect 
of the radial field. 
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Figure 1. Radial mode structure of the fastest growing n = 
branch MRI modes with 9 = tt/2,z = 0, S = 2.0 • 10"^ and 
Sg = 10 (solid), 10^ (dashes), 10^ (dots/dashes). The (purely real) 



growth rates are "fm = 
km = 36.0,30.0,20.6. 



0.39, 0.31, 0.19 and associated wavenumber 



(7) 

(8) 



where 

^2 

p =—((7 +7 )(a +7 + SbIcf cos 6), 
-til 

2_ ^ r i , 3 I 2, . N 

V =— — (7 +037 + a27 +ai7 + aio), 

and 

E\ — (cr^ + 7^)(cr^ + 7^ + Sb'JO-^), 

0,3 = Sbct^Icos^ 6(cot^ 6 — I) + sin^ 0], 

a2 = 2[cr^(l + cot^6l) + fc^] +2i5B(7^sin26lfc(l - cot^f), 

ai = ^flo-* [cos^ 61 (3 - 2fcV"^) + 1 - cot^ 9] - Sicr^fccot 9, 

ao = cr^[cr^(l - 2cot^ 61) - 2fc^]. 

Here a = kB sin 9 is the vertical Alfven frequency and Sb = 
Zv^/B^ ( = 3/^Bno/(B^7imi) in dimensional form) ~ fi/vu 
is a measure of the relative sizes of the anisotropic pressure 
force and Lorenz forc^ll- Here n is the number density, rrii 
the ion mass and /3 the plasma-beta. 

To obtain solutions to equation ([S}, we choose our 
boundary conditions to be an impenetrable wall at ro so 
5ur( ro) = and the requi rement r^Sur decays at infinity 
(see iFurukawa et al.l (|2007l ) for an inviscid treatment that 
includes the shear singularity at the origin). 

In this case, the eigenfunctions of equation ((6| are mod- 
ified Bessel functions of the second kind Ki„ (pr) with argu- 
ment pr and order v (|Watsonlll944l : lAbramowitz fc StegunI 
1 19641 ). The spectrum of solutions is discrete and we index 

V with n, Vn- In the special case when p (or 7 from equa- 
tion (O) is real (in general it is complex) the problem is 
Sturm-Liouville and Vn is an infinite ordered set of eigenval- 
ues vq < vi < V2 ■ ■ ■ < Voo ■ To determine «n and therefore 



* Us ing a Landau fluid clo s ure and including the effect of colli- 
sions iQuataert et al ] liool); ISharma et al. I liooi) were the first 
to show the dependance of the MRI on f2, vu and /3. Because the 
closure for the pressure anisotropy differs between the collision- 
less and Braginskii regimes, the transition between the pressure 
anisotropy driven MRI, and the MHD MRI scales, in dimensional 
form, as (yal^^fi^l'^ there, and {yiilSX)fi here. 



7, it is necessary to solve 

K™„(P) = 0. (9) 

From solutions to this and equations (O and ([SI, the full 
set of fiow, magnetic and pressure fields can be constructed. 
Appendix [A] In Fig. [T] we show the radial structure of &Ur 
when 9 = 7r/2. 

In general, determining 7 must be done numerically. 
However, for the most physically relevant magnetic field 
configurations, the problem becomes, in part, analytically 
tractable. These purely azimuthal field = (in 

the galactic plane), a purely vertical field 9 = -k /2 (a super- 
galactic field) and a slightly pitched field 9 (very slightly 
out of the galactic plane), exhibit categorically different be- 
haviour such that — >■ 0, 7r/2 are singular limits. 

In the first case, linear perturbations are damped; in 
the second, the system is unstable to the MRI but the 
Braginskii viscosity reduces the maximum growth rate 7m 
below the Oort-A value maximum |dlnf2/dlnr/2| = 1/2, 
(|Balbus fc Hawlevl[l993); in the third, the system is unsta- 
ble with 7m — » ^/|dln^^/ dlnr| = -\/2, even for asymptoti- 
cally small 9 (|Balbusll2004l ). We present the details of these 
calculations now. 



3 AZIMUTHAL FIELD, 61 = 

The stability of inviscid axisymmetric perturbations to a 
purely azimuthal magnetic field in the presence of a shear 
flow are well known. When q < 2 and the magnetic field 
B = Br~''e^ has d > —1 the system is always stable 
(as ours is). When only one criterion is met, depending on 
the form of the fields, the system may still rema i n stable 
by the modified Rayleigh cr iterion (|Ravleighlll916l : iMichaell 
I1954I : IChandrasekhail Il96lh . (It is worth noting however 
that global no n-axisymmetric perturbations are unstable 
(|Ogilvie fc Pringlc. 1996).) 

Do these results for our inviscidly stable system, TZ{'y) = 
0, persist in the presence of the Braginskii viscosity? The 
answer is no. 

Setting 9 = in equations (|6} to ([SI and making a 
change of variables, 5ur — >■ r~^^^Sur, we obtain the simple 
expression 



7 



2 d^SUr 



where 

3 1 



(7^Qi + 7Q2 + Q3) 5ur = 0, 



(10) 



Q 



4 



+ k\ O2 = Sb 



B^k^ 



2k^ 



(l + B^ 



We multiply equation (|10|l by the complex conjugate of Sur, 
5ur^ and integrate between ro and infinity. Boundary terms 
vanish, so 



dSur 
dr 



+ Ql\SUr 



+ 'yQ2\SUr\^+Q3\SUr\^ =0, (11) 



where |5urP — dr SurSur^ ^ is non-negative, as is 
(dSur/dr)'^. 

Equation (|ll|l is a quadratic in 7 whose roots de- 
pend crucially on Sb- When Sb = (the inviscid limit), 
7 is purely imaginary (neutrally stable, travelling waves), 
whereas when Sb > the result is quite different. In this 
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Figure 2. Left panel: Growth rate of the n=0 (fastest) branch of the inviscid global MRI as a function of the vertical Alfven frequency 
kBsinO for several values of B (time and length are non-dimensionalised with respect to ro and f2o). In the weak field limit, we recover 
the behaviour of the local MRI (dotted line). For all figures, the top-to-bottom order of the varied parameter corresponds to the top-to- 
bottom order of the main-panel curves. Right panel: Maximum growth rate 7m as a function of Bs'mO. As B 0, 7m asymptotes to 
the local Oort-A maximum = 1/2 (dotted line). Inset: The wavenumber km at which 7m occurs as a function of B along with the local 
maximum \/3/i (dotted line). Note: In the inviscid weak-field limit, two configurations with the same vertical magnetic field will exhibit 
identical behaviour - Section 15.31 



case TZ{'y) < V Sb and the only question is whether per- 
turbations are purely damped, or damped and travelling. 

It follows that if the system is stable in the absence of 
the Braginskii viscosity, it remains so in its presence. 



4 VERTICAL FIELD, 9 = 
4.1 Inviscid MRI, Sb =0 



: 7r/2 



For simplicity, we start by considering the inviscid limit of 
a purely vertical field B = B^e-, i.e. 6 = tt/2,Sb = 0. In 
this case equations ((6]) to ((8| simplifjjf]: 

Ki„„(fc) = 0, (12) 

with 



2k^ 



+ 1 



(13) 



If Dn < 0, Kii,„(fc) has no nontrivial zeros and it follows 
that for an instability > (hence our sign convention in 
equation (O). This implies 7^ is bounded from above by a^. 

Numerical solutions to equations (|12|) and (|13|l are ob- 
tained using Newton's method (implemented in Mathemat- 
ica). There is an unstable solution (and three stable ones) 
which is shown in Fig. (2] We find 7 is real, positive, and of 
order the shear rate. The instability is the global MRI, and 
the n — 0,vo branch has the largest growth rate. 

The fastest growing modes have k ~ 1/B, so in the 
strong field regime B < l,k > 1 the mode is large scale 
(physically, smaller scale modes are suppressed by magnetic 
tension) and is peaked away from the inner boundary. The 



^ Equations II12II and l|13| l correspond to equations (3.4) and (3.5) 
of ICurrv et al.l b994\ with a = 1 and equations (30) and (31) of 
iDubrulle fc KnoblochI \l99± who have already solved this prob- 
lem. 



result is that 7m is reduced below the Oort-A maximum that 
occurs when i? ~ 1 /k <^ 1 and the mode is localised at r = 1 
ICurrv et al.lll994l ). 

In this weak field regime the problem is amen a ble to 
asy mptotic analysis. It has been sh own bv lCochranI (|l965l ) 
and iFerreira fc Sesmal (|l97(]| . \200i ) that in this limit, the 
zeros of Ki^^ (fc) are given by 

«„~fc + s„2"'fc4 + ..., (14) 

where Sn = an23 , an is the m odulus of the n*** real negative 
zero of the Airy function Ai jAbramowitz fc StegunI 11964 ) 
and omitted terms are of the form k'' with 6 < 1/3. This 
result can be used to find the asymptotic-in-fc form of 7 by 
inverting equation (|13|l to form a bi-quadratic 



7^ + 2 



2 I 
cr + 



2 I 2/2 
7 + (T I (J 



into which we substitute equation (|14[) . Retaining the first 
two terms in Vn, we solve exactly for y'^ 



7 



+ 1 



Snk 



-1 ± 



l + 4a2 



(16) 



— Snk 3 

in which the positive root corresponds to the instability. 
Numerically we find the fastest growing mode (for a given 
B, Sb), 7m and the wavenumber at which it occurs fcm, obey 
dym/dB,d{kmB)/dB <0. 

To draw an analogy with the local approximation (see 
Section I6.3|l . n that indexes the number of zeros in the 
domain is like a radial wavenumber. For large enough n, 
the solutions are plane waves. To see this, we apply a 
WKB analysis to equ ation ^ using the small p a rameter 
kr/v n = x/vn 1 IDubrulle fc KnoblochI 1 19931 : lOgilvid 
Il998l ). We have 

___ + —SUr{x) = 0. (17) 

Demanding that Sur is real, this has solutions 

5itr = %/scos(i)n Inx), (18) 
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and to ensure cos(i;nlna;) = at r = 1, it must satisfy the 
boundary condition 



(19) 



Wn In fc = I n + - I TT, 

which determines the spectrum of solutions. 

Combining equations p3|l and (|19|l we have 



7 



(20) 



where a — k^{lnk)^ /n^-rr^ and we have neglected factors of 
1/2 « n. 

In the large n limit, on small enough scales, the mode 
structure given by equation (|18|l is especially simple. Revert- 
ing to r as our radial coordinate, we expand equation (|18p 
about r = ri ~ 0{1). To leading order we find 



Sur^A{l + ~j cos + ^ 



(21) 



where A = v^fcrT and ^ = Wn In kri. 

These solutions describes rapidly oscillating modes with 
frequency Vn and a slowly varying amplitude cx A. For suffi- 
ciently large Vn, the solutions are plane waves whose growth 
rate decreases with n. 



4.2 Viscous MRI, & / 

Allowing for viscosity whilst retaining a vertical field, equa- 
tions (0 and (ISJ again reduce to a simple form: 



2 I 2 

7 -I- o- 



2fc2 



(7^ - 



7^ + -f Sb'JO'^ (7^ + o"^) 



-HI 



(22) 
(23) 



Unlike the inviscid case, we can no longer guarantee the re- 
ality of p as 7 can, and indeed sometimes does, take complex 
values. The problem is generally not of Sturm-Liouville form 
and so the roots of Ki„„ (p) are complex. 

Numerically we again find four branches of which three 
are stable, and one unstable. The unstable mode (the only 
one of interest) is real, so if 5R(7) > then 7 G 5R. In this 
case, the problem is Sturm-Liouville. The unstable branch 
can be traced from the inviscid MRI (and is maximized for 
n = 0,vo), Figs. [H [3] and [5] 

Asympto tic e xpressions for 7 and k can be found using 
the results of ICoch ran (196^ and lFerreira fc Sesnia (|2008i ') 
for the (complex) roots of equation 



Vn ~ Tt/ (In 2 - 7Euler ' 
Vn ~ p + an2'ipi + 



Inp), 



P G 
P G 



bl « 1, 
\p\ » 1, 



(24) 
(25) 



where 7Euier — 0.58 is the Euler constant. 

For Sb S> 1 we combine equations (|22|l and (|23|) to get 



fc4 k^ 



Re ^7 - 



1 + v;- 



+ 2B' 



+ Re" 7 



1 + 



(26) 



where Re = SbB is the inverse of the Reynolds number. 



We expand equation (|26p in 7^/A:^ ~ e <C; 1 and solve 7 
70 + £71 + . ■ . order by order. To lowest order we find 



70 



Re 



1 + vl 



B' 



and to 0{e 
1 



71 = 



£2 



and so 



7~Re 



l + vl 



70 + 2Re 



--B'\- 



1 + vl 



+ B' 



k^B-2 



70 + 2Re 



■ + B- 



(27) 



Numerically we find fcm ^ 1 which, combined with 
equations (|22p and (|24p . implies Vn — — vr/lufc. It follows 
that to lowest order 



7n. = Re (2 - B) 



5i3 > 1, (28) 



and, differentiating equation H27|) with respect to k and ne- 
glecting leading order logarithmic variations, 



B' 



Btv 



R2 

2 -f — Re (In Re" 



■1x3/2 



S'b>1. (29) 



For Sb,B 1, equation (|25p applies and so 7 is gov- 
erned by 



(7' +aY+[2^ k-^ {SB-,a^ - 2)] (7' - a') = 0, 

where we have used Vn — k ^ 1. We expand in 1/fc^ ^ 1 
so 7 = 70 + fc^^7i + . . ., and again solve order by order. To 
lowest order we find 70 satisfies the inviscid equation 



70 = V -(cr2 + 1) + v/l + 4cr2, 
and at C»(/s"^) 

-7o (1 + SBloa^) +a^l- 5'fl7oo-2) 



71 



4(70 + o-)'^ + 70 



Sb-^I, (30) 



, Si3 < 1. (31) 



To see the effect of variations in n on 7 we perform 
a WKB analysis in which we order Wn ^ 2> 1 ^ Sb so 
\np/vn is a small parameter - see Section r4.1l The boundary 
condition for the viscous modes is then 



Vn In 



k 1 



ISbJo^ 
2 72 -h o-s 



Combining this with equation (|23|) . we expand in Ss <^ 1 
so 7 = 70 -I- 5'b7i 4- . . .. To lowest order 70 is given by the 
inviscid equation (|20|l . and to 0{Sb) 



-7ocr^V7^ + ^ ^ n 

71 = < 0, 

4(70 + o-)3 -I- 2q7o 



«n>l,S's<l. (32) 



In the various limits, our asymptotics confirm the nu- 
merical results that 7 G 5R, and O^/OSb < 0. Contrasting 
the results of this section with the inviscid case when B ^ 1, 
we have 



Sb 
Sb 



0, 

■ 00, 



7m 
7m 



1/2, 
Sb'B- 



B''\/374, 

Sb^B-^IuSbB^^^. 
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Figure 3. Growth rate of the n=0 (fastest) branch of the global MRI with Braginskii viscosity for a two different weak, B = 10"'^, 
field configurations and a range of Sb- Left panel: Vertical Field: For fixed k, 7 decreases as Sb increases and is bounded above by the 
inviscid growth rate (dotted line), but the critical wavenumber above which 7 < is independent of Sb- Right panel: Slightly pitched 
field; In contrast to 8 = 7r/2, for a given k, the viscous case jISb 7^ 0) > 7(5b = 0) inviscid case (dotted), but is less than the local 
Lorentz-force free, or MVI, limit (dashed). As Sb ~^ 00, km — > but, unlike the vertical case, 7m — > \/2- 



5 PITCHED FIELD, e < 1 
5.1 Ordering assumptions 

When the magnetic field has both a vertical and an az- 
imuthal component, the perturbed Braginskii stress tensor 
exerts an azimuthal 'tension' force on separating plasma ele- 
ments. For arbitrary 9, the system is neither Sturm-Liouville 
(equation ((6| is complex) nor is it amenable to the kind of 
polynomial inversion used in Section |4l 

However, assuming 6 <^ 1 matters simplify consider- 
ably. Physically, this choice of pitch angle represents the 
most r ealistic non-isolat ed galactic magnetic field configu- 
ration (|Beck et al.lll996l '). Mathematically, as we now show, 
it is a singular limit that constitutes the stability threshold 
between actively damped modes, Section O and an unsta- 
ble configurations that grows even faster than the inviscid 
global MRI (to which the following stability threshold also 
applies) 

Damped -.6 = 0, Unstable : 9 ^ e. 

Here e ^ 1 is a small parameter with respect to which we or- 
der the pitch angle and the remaining quantities in equation 
©: 7,fc,d/dr,B and Sb- 

Letting 9 ~ e we can retain only the first few terms in 
a series expansion of our trigonometric functions 

cose = 1 - y +0(6i^), sm9 = 9-^ + 0{9''). 

To include their effects, we assume a strong magnetic 
field B ~ 0(1) and Braginsku viscosity Sb ~ 0(1). To 
retain vertical magnetic tensions a — kB sin 9 to order 1 
(failing to do so removes the high wavenumber cutoff and 
leaves the equations ill-posed), we set k ~ 1/e. Then, in 
anticipation of unstable modes that grow at the shear rate, 
we order 7 ~ 1 too. Balancing terms in equation ([SJ, we find 
d/dr ~ p ~ fc ~ 1/e. 

In summary, our orderings are 

9^e, 7~o-~B~S's~l, fc~-i~e"\ (33) 

dr 



We apply these scalings to equations ([61) to ([S]) and re- 
tain the lowest order e terms. We find Sur is still governed 
by Bessel's equation and, since p = fc, 7 is determined by 
equation (|12|) and 

i;^ = --^-(7'* -I- 637^ -I- 627^ + fei7 + feo), (34) 
where 

El = (a^ + 'y'^){o' + 7^ + Sb^o^), 
h = SB(y^9^'^ , 

62 = 2fc^ (1 + B^) ~ 4iSscr^A;0"\ 
bi = -SB(7^k'^ (2 + B^) - Sicrfc^B, 
60 = -2cr^fc^ (1 + B^) . 

We start by solving the inviscid problem. 

5.2 Pitched inviscid MRI. Sb =0 

In the inviscid weak-field regime, we recover the vertical 
weak-field instability of Section l4Tl and Fig. [21 B appears via 
the vertical Alfven frequency only so fcm(6' — 7r/2)/fcm(e ~ 
tilted) = 1/0 S> 1 so the unstable mode will be confined to 
a boundary layer of width ~ 1/{B9) ~ 0(e^)- 

If B is small but finite (implicitly all orderings are sub- 
sidiary to equation d33}) the governing dispersion relation is 
the complex quartiqj 

(7^ + o-'^ f + 2(1 + B^)(7^ - a^) - 8iB(J7 = 0. 
Writing 7 as a series in B ^ 1 we find 



^ Solutions to this p olynomial are considered in detail in 
Curry fc Pudrit3 l|l995h and further discussion can be found in 
KnoblochI l ll992l'l . 
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Figure 4. 5R(7m), Q(7m) (if it exists), and the corresponding vertical Alfvon frequency fcm-BsinS for two field configurations. Left panel: 
Vertical field: 7m S 5R and kmBsind vs. Re~^= SbB^ for a range of B. These are bounded above by the B <g; 1 local limits of 1/2 and 
-y/3/4 (dotted lines), and well matched by the asymptotic results (dashed lines) given by equations II28I I and II29I I. Right Panel: Slightly 
pitched, e = 10-3 : 5R(7m) and kmBsiYiO are bounded from above by the local -B <§C 1 limit of \/2 (main panel, top dotted line) and the 
inviscid limit ^3/4 (top inset, dotted line). The imaginary part of the growth rate 3(7m) (bottom inset) is well described in the small 
and large Sb limits by equations l|36| l and l|48|l . Note that d^jmiB ~ l)/dSB < (bottom curve, main panel) demonstrates the combined 
effects of hoop tension and the Braginskii viscosity. 



with a maximum value 
' 2 V 4 

which occurs at km = \/3/4:B~^ , Fig. U) 



(36) 



5.3 Pitched MRI with Braginskii viscosity, Sb 

In the presence of the Braginskii viscosity we find a com- 
plex instability whose real part exceeds that of its invis- 
cid counterpart. Numerical solutions described by equa- 
tions ([12} and (jSlll are shown in Figs. |3l |4] and [5] When 
-B — !> 0, Sb oo the growth rate tends to the shear rate 
5R(7) — > \/2d In Q /d In r = \/2 and, as before, the fastest 
growing modes occur for n = 0. 

Asymptotic expressions for 7 and k can be found by 
substituting equation (|14p into equation (|34|l . For B ^ 1, 7 
is governed by 

{-/^ + a^f + 2(7^ - a^) + 5i3cr^7(7^ + a" - 2) = 0. 

Differentiating this with respect to cr {k, B and 6 appear 
together only in one combination and so a is treated as a 
single independent variable), and setting A'y/da = 0, the 
stationary points of 7 obey: 



SbY + 27 + 2Sb{(t - 1)7 + 2(<t' - 1) 



0. 



(37) 



By considering the branch a^n that maximizes the stationary 
value of 7, and introducing C,-m = 

Ss7m(2-7m)+2(l-7m) 



(38) 



2(7mSs + 1) 

and substituting this into equation H37p we obtain a poly- 
nomial whose solutions describe 7m: 

7m-47m+4(l-4SB')7,i+45s'(27m+5s'-57^L) = 0. (39) 

Because it is a sixth order polynomial, there is no gen- 
eral formulae for its roots. However, the presence of the 



asymptotic parameter Sb recommends a series solution. The 
power of the expansion parameter Sb depends on whether 
it is large or small. 

Taking Sb ^ 1 first, we write 

7m = 7m,o + 5b7™,i + C('S'I), (40) 

Cm=Cm,0 + SsCm,l+O(S|). (41) 

Substituting equation into equation and solv- 
ing order by order, we find, to lowest order 7„i,o = 1/2. To 
next order 7m, i = 3/32 and so 



^-=2(1 + ^-16 



Si3 < 1. (42) 



Substituting this into equations (|38|) and (|41|l , we find cTm = 
•\/3/4 (1 — Ss/48), so the most unstable mode is 



km — 



Sb 



1 



Be\i A 

Now taking Sb 3> 1, we write 



Sb<^1. (43) 



7m = 7m, + 'S'g^''^7m,l + 5a^7m,2 + 0{Sg^^'^), 
Cm = Cm,0 + Sg^^^(^m,l + S b^ Cin,2 + 0{S^ 



3/2n 



(44) 
(45) 



Following the same procedure, to lowest order 7m, = 
±\/2, and we take the positive root corresponding to the 
instability. At first order we obtain no information, but at 
second order we find that 7m, 1 = ±2^^* . Taking the 
negative root (see Fig. [S| we have 



7m 



^/2 1 



9 



Sb»1. (46) 



Substituting this and equation (|45|) into equation (|38|l . to 
lowest order ^m.o = 0, and at next order Cm,i = 2^^^ /3^^'^ so 
am = (2V9)^/»S-^/^and 



1/8 

0-1/4 



Si3 > 1. (47) 
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Figure 5. Direct comparison of the 7ni (main panel) and 
fcmSsinS (insets) behaviour of the unstable mode for d = 10~^ 
(top solid line and inset) and 6 = it/2 (bottom solid line and inset) 
with Sb for a weak mag netic field, B = 10"^. The dotted lines 
are the local maxima, and the dashed lines are the asymptotic 
scalings given by equations I I28II and I I46II for -y^, and equations 
|[29t and gTjl for fcmBsin0. 



Now considering the effect of finite magnetic fields, 
5R(7m) (and 9(7m) which is now present) increase with Ss- 
Guided by the numerical results, we assume I7I ~ 1 and take 
the leading order balance of terms oc Sb — >■ 00 in equation 
We find 

7^ (cr^ + B^) - AiaB-f + a^{a'^ - - 2) = 0, 
and the resultant growth rate is 



7 ■ 



2iBa + a ^ (B'^ - a^) - 2{B - a) 
52 +cr2 



, 5 < 1, > 1, (48) 



which agrees well with the numerics, Fig. U) 

To see the effects of n on 7 we adopt the same WKB 
approach and orderings as Section [4. 21 where Vn is now given 
by the appropriate limit of equation (|34|) . We assume B <^ 1 
and expand in 5s ^ 1 so 7 = 70 + SbJi + . . .. To lowest 
order 70 is given by equation (|20p . and to 0{Sb) 



^ joa (2a - 7() - ^ ) ^ 
4(70 + cr)^ + 2a7o 



(49) 



so the Braginskii viscosity increases the growth rate beyond 
the inviscid limit. 

Contrasting the results of this section with the inviscid 
case, for B ^ 1 we have 



Sb 
Sb 



0, 



K(7n 

K(7„ 



1/2, 



iOB)-'^, 

{Bes'J')-\ 



with the caveat that fcm — >■ only in the formal limit B — >■ 
too (see Fig. |4]for the effect of finite B at large Sb)- 

Considering the effects of hoop tension, we find ^(7) 
is an order of magnitude greater in the viscous limit. For 
example, we find numerically that ^{^{Sb ~ 10*)/7(S's = 
0)) ~ 40 for B = 10~^ (which agrees with equations H36p and 
(|48p to within half a percent for this and all other examples 
with B < 0.5, Sb > 1.) 



6 DISCUSSION 

6.1 Physical mechanism 

We have seen that the orientation of the magnetic field cat- 
egorically determines the behaviour of the system. How is 
this to be understood physically? 

The most informative explanation comes in the weak 
field regime where the role of the magnetic field is twofold. 
Firstly it facilitates the generation of pressure anisotropics 
proportional to its rate of change, or equivalently (in the col- 
lisional regime) oc Sbb : Vu. Because collisions are involved 
this is a necessarily dissipative process. Equations ((T)) to 
yield an energy conservation law 



d_ 



+ 



3i^B {\bb : V-uf 



Sub 



/dln_B 
V dt 



where < ■ > are volume averages. 

That is, the Braginskii viscosity damps any motions 
that change the magnetic field strength. This is the first 
role of the magnetic field. 

The second role of the magnetic field is a geometric 
one. Assuming anisotropics do arise (from changes in the 
magnetic field strength), the field's orientation dictates the 
projection of the anisotropic stress onto the fluid elements of 
the plasma, thereby affecting their dynamics. As we now ex- 
plain, this fact is crucial in determining the stability, or lack 
thereof, and the role global effects have on a differentially 
rotating magnetized plasma. 

When the field lines have both vertical and azimuthal 
components (it is pitched) fiuid elements at different heights 
can exert an azimuthal stress on each other. The sign of 
this stress can be either positive or negative. If the mag- 
netic field is unstable, so its rate of change is positive, the 
anisotropic stress acts to oppose any azimuthal separation of 
the two elements (this can be identified as the fiuid version 
of the stress responsible for the microscopic mirror force). 
In this case, like the MRI, velocity perturbations to fiuid el- 
ements at different heights increase (decrease) their angular 
momentum causing them to move to larger (smaller) radii. 
In a system with an outwardly decreasing angular velocity 
profile, this leads to an azimuthal separation of the fiuid el- 
ements. The associated magnetic field growth ensures the 
stress is of the right sign to oppose this separation and this 
transfers angular momentum between them in a way that 
facilitates furth e r (rad ial) separation; i.e. an instability (see 
iQuataert et al.l (|2002l ) for a physical explanation including 
a spring analogy). This is the second role of the magnetic 
field. 

In conjunction, the two roles explain our results. In iso- 
lated field geometries where there is no projection of the 
stress onto fiuid elements at different heights, the Bragin- 
skii viscosity does not give rise to an instability ; its only ef- 
fect i s dissipative (|Kulsrudll200"5l : lLvutikov|[2007l : iKunz et al.l 
l201lf ). This accounts for the damping of perturbations in 
the azimuthal configuration. It also accounts for the re- 
duced growth rate of the vertical MRI (that depends on 
the Maxwell, not the Braginskii, stress), along with the in- 
creased radial exten t of the mode aw ay from the region 
of maximum shear (|Currv et al.lll994l ). When the field is 
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pitched its dissipative effects persist, but the free energy 
contribution from the differential shear flow will always be 
greater, leading to an instability (the dissipative effects can- 
not change the stability boundary, just the growth rate). 
Of course the total, equilibrium plus perturbed, field energy 
will decrease if viscosity is present (and viscous heating is 
not) but, from the perspective of generating turbulence or 
transporting angular momentum, this is a secondary conse- 
quence. 

Now considering finite magnetic field strengths (and 
therefore hoop tension and viscous curvature stresses) we 
find these modes become over-stable. The variation of the 
travelling wave component of the unstable mode is simply a 
combination of these two effects, to varying degrees. 

6.2 Astrophysical example 

The mathematical results in this paper should be appli- 
cable to any collisional magnetized disc with a decreas- 
ing angular velocity profile. However, it is helpful to pro- 
vide a physical example where our analysis holds. Be- 
cause our theory is developed for a disc with an angu- 
lar velocity profile oc we choose the ISM of a spi- 
ral galaxy where a range of analy t ic and numerical stud- 
ies a l ready exist jKim et al.l 2003| : Kitchatinov fc RiidigeJ 



consequently the mean free path is 



Wang fc Abej|2009l ). 



20041: iDziourkevitch et all |2004| : IPiontek fc QstrikeJ bool 



Further motivation t o study this system comes from 
ISellwood fc BalbusI j 19991 ) who have argued that in the HI 
region outside the optical disc of a spiral galaxy, the velocity 
dispersion measurements of ~ 5 — 7 km s~^ may be driven 
by the differential rotation of the disc, mediated by the MRI. 

On this basis, we consider the warm ISM in the qui- 
escent regions of a typical spiral galaxy where the plasma 
is generally subject to a weak magnetic field and an out- 
wardly decreasing angular velocity profile. It is also mag- 
netized, and so subject to the effects of the Braginskii vis- 
cosity. This last feature can be seen, and the other rele- 
vant parameters estimated, by a dopting the following set of 
fiducial parameters for the IS M (jBinnev fc Tremainelll98"^ : 
iBeck et al.ll 19961 : iFerrierdbOOlh . These are in agreement with 
NGC1058, the well s tudied face-on disc galaxy considered by 
ISellwood fc BalbusI (| 19991 ). 

Reverting to dimensionalised units for clarity, we have: 

• Particle number density (ion and electron are the same) 
n ~ 0.3 cm~^. 

• Temperature (we assume ions and electrons are in ther- 
mal equilibrium) 

r, ~ 5 X 10** K; 

consequently the ion thermal speed is 



Vthi 



,= l^] ~ 3 X 10" cm s-^ 



1/2 



(Ti is in erg.) 

• The ion-ion collision frequency (in seconds, assuming n 
in cm~^ , Ti in K and the Coulomb logarithm In A = 25) ifl 



1.5nT, 



-3/2 



' 4 X 10" 



Vthi 



7 X lO" cm. 



The typical rotation rate of a spiral galaxy is 



n ~ 5 X 10"^" s"\ 

A typical value for the outer edge of the optical disc where 
the turbulence cannot be generated by stellar processes is 

r ~ 3 X lO'^^ cm, 

and even within the optical disc, at the corotation radius in 
between the spiral arms, magnetized shear instabilities may 
be important (A. Shukurov - private communication). 

Outside the optical disc a reasonable value for the vertical 
scale-height is 

H ~ 10^* cm, 

and so the disc is thin. The measured system-scale rotation 
(not turbulent) velocity is 

[/ ~ 2 X 10^ cm s~\ 

• The observed mean magnetic field strengths vary be- 
tween galaxie s but on the lowe r end of the scale are, at the 
present time (|Beck et al.lll996l ) 



B ~ 8 X 10^ cms"^ 



present. 



However, if the MRI is the dominant turbulence gener- 
ating mechanism in the ISM, this value must represent 
the saturated state of the magnetic field. Assuming, as 
one must, that present field strengths have been ampli- 
fied over time, at some e arlier ti me they were weaker, 
e.g. iMalvshkin fc KulsrudI (|2002l ): iKitchatinov fc Rudigerl 
(|2004l V 

To ensure the most unstable modes exist at scales 
> Amfp so our theory is fluid-like we adopt the follow- 
ing value for the historical 'initial' field str ength and lay 
aside the problem of where this ca me from (|Kulsrudlll999i : 
[Brandenburg fc Subramanianll2005l ). 



B . 



initial. 



If we considering a plasma in this era, or the ISM in a galaxy 
where the magnetic field is not saturated, the plasma beta 
is 



/? = ~ 1.3 X 10^ 



the ion cyclotron frequency is 



The full expression for the ion-ion collision frequency (ion- 



U-jvneB 6 -1 

iUi = 1 / ~ 2 X 10 s ; 

V riii c 

(e is the elementary charge, c the speed of light) and the ion 
Larmor radius is 

pi = ~ 1.5 X 10 cm; 

and so the plasma is magnetized uji 3> I'a (and will become 
even more so if the magnetic field becomes stronger). 



electron collisions are sub-dominant for a thermalised plasma) is 
given by ua = i^/wne* \nA/3ml^^T^^^ where e is the elementary 
charge and T is in erg jBraginskiill 19651) . 
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• The magnetic Prandtl number Pm is huge (and so we 
can take the Braginskii viscosity to be the only significant 
dissipative process 



where 



P„ = — ~ 7.5 X 10"''- 
V 



7.5 X 10' 



where rj is the coeiEcient of resistivity and T is in K. 

• The dimensionless parameter Sb ~ SuB^i/ {nrriiB^) ~ 
130, /vii is 

Sb ~ 1.1 X 10^ ~ 42, 

for the conditions above, so the Braginskii viscosity plays a 
role. 

Our model of a gravitationally constrained, differen- 
tially rotating disc made up of an isothermal, incompress- 
ible, magnetized plasma fluid is consistent. Neglecting struc- 
ture in z requires us to restrict our analysis to the galactic 
plane by considering vertical scales much less than the scale 
height of the disc ~ Cs/fi ~ where Cs ~ t^thi is the isother- 
mal sound speed. 

Formally, both this and our remaining model assump- 
tions can be expressed as a hierarchy of time-scales which are 
well satisfied for our set of parameters (we restrict attention 

to 7m), 



1 



1 1 ro 



In order of increasing periods of time, these scalings corre- 
spond to: the plasma being magnetized; the plasma being 
coUisional, i.e. a fluid; the model being uniform in the verti- 
cal direction; and the disc being thin, i.e. rotationally rather 
than pressure dominated. 

The parameters in this section describe a physically re- 
alistic regime in which our analysis is valid, and potentially 
important in explaining the gas motions in parts of the ISM. 



6.3 Relation to the local approximation 

Because of its widespread use, we briefly comment on the 
relation of our results to the local analysis. In this ap- 
proximation the coordinate system is Cartesian, the shear 
is modelled linearly, and the perturbed quantities are as- 
sumed to vary rapidly in the radial direction (with respect 
to the background variations) so they may be written as 
5u{r) = (5uexp[i(Zr -I- kz) -\- ')t], Ir 3> 1. Under this as- 
sumption, a WKB analysis ignores terms proportional to 
1/r and replaces d/dr by il. The loc al dispersion relation is 
(llslam fc Balbus|[2005l : lFerrardl2007l ): 



7" + di7^ + d27^ + 0(37 + d4 = 0, 



(50) 



* The coefficient of the Braginskii viscosity is given by i/b = 
0.96nmiv'^^^^/uii and 77 = 3y/2me /TTc'^ In Ae^T"^/^ where nie i s 
the electron mass and T is in erg llSpitzeJll963 : lBraginskiilll965^ . 



di =Sbo' 



2/^ -f fc^ cos^ 



dz =SBcr a 



P + k^ J ' 
2 + fc^ cos^ ( 



d4 —a^ i — 2 



P + k'^ 
k' 



P + k2 



P + k^ 



In the A: > /, ; > fc limits for 6 = 0, 7r/2 the local 
and global results are identical (modulo different boundary 
condition-dependent restrictions on the spectrum of allowed 
modes in the latter case). 

However for 6 ^ 1, B < 1 the analyses differ. Locally 
we find unstable modes have 7 G K, whilst globally 7 is 
complex. Whilst in the inviscid limit 9 (7) ~ B can often be 
neglected. In the highly viscous limit ^(7) ~ aB/{a^ + B^) 
and, as shown in Fig. [J] this can be over an order of magni- 
tude greater. This difference may prove important (from a 
modelling perspective) for viscous systems which could have 
been treated locally, were they inviscid. 

Furthermore, inconsistent with the global picture, for 
9 — 0,n/2 and I — the local analysis neglects the Bragin- 
skii viscosity. This can be understood from 

5{bb : Vu)g = 60 ■ (60 • ViSUrBr)) = — , 6* = 0, 

r 

5(bb : Vu)g ~ 62 • (bz ■ V(5uzez)) = ikSuz, 9 — 7r/2, 



whereas locally, assuming V 
have 

5{bb ■ Vu)l = 0, 



5u = i{k5uz + ISur) = we 



= 0,7r/2. 



(Here the subscripts G and L stand for global and local.) 
Physically, in the global component of the fiow is 

projected along the magnetic field direction by the curvi- 
linear geometry (azimuthal case) and the demands of the 
global incompressibility condition (vertical case), and so \B\ 
changes at linear order, thereby activating the Braginskii 
viscosity. This is not so in the local case. 

One final clarification is worth noting. The local de- 
scription predicts that the fastest growing unstable modes 
(9 7^ 0) occur as i — 0. Formally this is inconsistent with the 
assumption Z 2> 1/r 7^ that went into deriving equation 
(|50[l . The local solutions obtained under the WKB approxi- 
mation are not self-consistent and should be described by the 
type of global solution we have constructed here. However, 
having determined the global solutions, we can confirm the 
n = branch corresponds to 7in, and so the local analysis 
is, at least qualitatively, correct in this respect. 



7 CONCLUSION 

The nature of the ideal MRI has been well established but 
how non-ide al modifi cations affect it re mains an active topic 



of research (iBlaes fc Balbiis 19 94; B albus fc TerauernI 
Kunz fc BalbusI |2004|: iFerrarol l2007l : IPessah fc ChanI 



2001 



200s 



Devlen fc Pekunlull201of r Of interest to us, and the subject 
of this study, has been the global nature of the Braginskii- 
MHD MRI. 
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Considering an isothermal, magnetized, collisional disc, 
perhaps the early ISM, we have found that for an azimuthal 
magnetic field with an asymptotically small vertical compo- 
nent, a singular instability emerges - the magnetized MRI. 
The growth rate of this instability ~ Q depends on Sb, a 
dimensionless combination of the temperature, shear rate 
and Alfven speed, and it is up to a factor of 2y/2 faster 
than its unmagnetized counterpart. In the limit where Sb 
is large, we determined the asymptotic maximum growth 



-1/2 



rate "/m — \/2 (^1 — A'^"' ' j and corresponding wavenum- 

ber km — {BB)"^ S g^^* . As the field increases in strength 
the Larmor radius decreases dpi/dt oc — (so the plasma 
becomes more magnetized) and as the azimuthal compo- 
nent of B increases, purely growing modes mutate into over- 
stabilities whose imaginary part also depends on Sb- Both 
of these factors may be important for the questions posed in 
Section [T] and for the turbulence that will eventually arise. 
In the early stages of the instability, before fully developed 
turbulence has set it, the instability will generate non-linear 
pressure anisotropies. What exact effect these will have is 
an open question, and so we finish with some thoughts on 
the subject. 

One major uncertainty in the evolution of magne- 
tized accretion discs is the effect of pressure anisotropy 
driven micro- instabilities whose growth rates are gen- 
erally well in excess of the MRI e.g. mirroi[f|and fire- 
hose jSchekochihin et al.1 120051 ). Although various mod- 
els for treating them exist , a complete first-principles 
theory remains outstand i ng dSchekochihin fc CowlevI l2006l : 
ISharma et al]|2006l . [20071 : iRosin et aLlbOHl ). It is crucial for 
accretion theories to determine the fate of these instabilities 
fluctuations over long (transport) time-scales. This is be- 
cause they determine the Maxwell and Braginskii stresses 
that, in turn, dictate th e angular momentum transport 
IShakura fc Sunvaevll 19731 ). Even if the pressure anisotropy 
is pinned at the marginal value for the microscale instabili- 
ties in some self-regularizing way, there are further compli- 
cations that must be addressed. 

Specifically, the viscous stress generated by the Bragin- 
skii viscosity can heat the plasma - at a rate oc [pj_ — py)^ 
(|Kunz et al.ll201ll ). Spatial inhomogeneities in the growth 
rate of the magnetic field, as one would expect in a turbu- 
lent system, will lead to inhomogeneities in the local pressure 
anisotropy, both in magnitude and sign. In regions of de- 
creasing field strength the firehose would pin the anisotropy 
at \{p± — P\\)/p\ = 2//3 and in regions of increasing field 
strength the mirror would pin it at \{p± — P||)/pi = 
The implication of this is that heating in regions of increas- 
ing and decreasing magnetic field strength could differ by a 
factor of ~ four and would occur on the decorrelation scale 
of the turbule nt cascade's viscous cut-off - w here the shear 
is max;imized (jSchekochihin fc CowlevI [200^ ). If differential 
heating of this nature does occur then one might expect 

^ In regions where dln_B/dt > slow- wave polarised modes 
(like tile MRI) are unstable to the mirror witii 7m ~ LJi{{pj_ — 
P||)/p)^. For the set of parameters listed in Section 16.21 this is 
^. a faetor 10^ slower than the MRI. However in 



2 ■ 10" 



many other contexts, e.g. the intracluster medium, the mirror 
can be up to a fact or ~ 10* faster than the macroscopic shear 
rate faeUingedbOOTl'l . 



the non-linear dynamics to be further complicated by mag- 
netized (and unmagnetized) temperature gradient instabil- 
ities, and the temperature dependanc e of both vb and the 
micro-instability thresho lds ((Balbus ,200ll : lQuataer"3 12OO8I : 
ISchekochihin et al.llioiol '). 

Understanding the rich interplay between these realisa- 
tions of magnetized plasma phenomena constitutes an im- 
portant and, as of yet, unsolved issue in astrophysics; specif- 
ically in accretion discs, but also in galaxies and galaxy clus- 
ters. The overall picture is a deeply interconnected one and 
the types of processes outlined above are probably pertinent, 
to some degree, to most magnetized astrophysical settings. 

To address these issues there is a need for both non- 
linear simulations that include a wide range of magne- 
tized physics, and a more complete theory of the trans- 
port effects of micro- instabilities. For the first task at least, 
this work should be usefu l for benchmarking global codes 
ISkinner fc Ostrikeij[2010l ). 
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APPENDIX A: LINEAR STABILITY ANALYSIS 

The linearised equations ([1} to Q are, in component form, 

(Al) 

(A2) 

(A3) 
(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

where Z = SsB^bb : Vn. Equations (|A1[) - (|A8[) form a closed set which we can combine into a single differential equation. 
Eliminating the perturbed magnetic fields, stress, pressure and the vertical component of the velocity field yields two coupled 
ordinary differential equations for 5ur and Su^ 

Co — —^+Cl5u^ = C2 — rT~C3 — :r~~+C!A5ur 



r 


dm 

dr 


+ io5Br 


B cos eSB^ dSZ cos^e 

— 2 — — ■ — — 3 oZ, 

r dr r 


^ SUr 

^dud, H 

r 


BcosOSBr ,sin26'^^ 

= h iooBa, + 3ifc — - — oZ, 

r 2 




= -ikSU 


+ iaSBz 


-ifc(f-3sin^ e)5Z, 


"fSBr 


= iaSur, 






1SB4, 


= iaSu^ - 


SBr 


B cos 9SUr 


r 


r 


-ySBz 


= iaSuz, 








1 d{rSur) 
r dr ' 




SZ 


3 


^cos^ 6 


fcsin20\ Sur 2 ,sin2S 

— 1 hifcsm Oouz + ik ous 

7 2 y r 2 



dr 

where 
Co = 



DoSua, = Di 



dr^ 

ASUr 

dr 



k^2 I f -i-Bcos( 

7 



dr 



ikSBcr^ cos^ 6 cot 9, 



(A9) 
(AlO) 

(All) 



-ifcSscr — - — , 

'T^Ht^ + f^^) — Sbct'^ sin^ 6, 

-{j^ + a') + Sbct^{Eo~ cot'' 



-{'•f+a) {k' + 



B cos e Lb cos e - i- 



Si3(7^_Bo(cot^6»- 1) 



(A12) 
(A13) 
(A14) 

(A15) 

(A16) 
(A17) 

(A18) 

(A19) 

Equations (|A9P and (|A10|l can then readily be combined into equation ([6]). Solutions to this equation are modified Bessel 
functions and so the functional forms of the perturbed fields are: 



Ci = 

C2 = 

C3 = 

Do = 

Di = 

D2 = 



7^ + (T^ + 'ySso'^ cos^ 0, 
—ik~^ Sb^o-'^ cos 6 sin 6, 



a I 2iB cos 6) + - - - 



2k'- 



+ tSb'Jo-^ cot 6E0 
k 



Eo = — ifc7 ^ sin 6 cos 9 + cos^ 6 — sin^ 



Sur — Ki„jj (pr) exp[ifc2 + ^/t], 

Di dSUr ^ D2 , 

Do dr Do 



5Br 

SB^ 
SB, 



i 1 d(r5ur 



k r dr 

l — OUr, 

7 

. a Di dSur 

7 Do dr 
a 1 d{rSur) 
^ kr dr ' 

(7^ + cr^) 1 d(r5w,.) 
k^-y r dr 



. a D2 ^ 1 I B cos ( 
7 Do r 



- SbB^ - sin 



7 



1 



5Ur 



£2 

Do 



d5ltr 

dr 



£0 . , sin 29 Di 

h IK =— 

r 2 Do 



SUr 



(A20) 
(A21) 

(A22) 
(A23) 

(A24) 

(A25) 

(A26) 



